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A comprehensive analytical theory for nonlinear excitations related to horizontal (longitudinal, 
acoustic mode) as well as vertical (transverse, optical mode) motion of charged dust grains in a dust 
crystal is presented. Different types of localized excitations, similar to those well known in solid 
state physics, are reviewed and conditions for their occurrence and characteristics in dusty plasma 
crystals are discussed. By employing a continuum approximation (i.e. assuming a long variation 
scale, with respect to the inter-particle distance) a dust crystal is shown to support nonlinear kink- 
shaped supersonic solitary excitations, associated with longitudinal dust grain displacement, as well 
as modulated envelope localized modes associated with either longitudinal or transverse oscillations. 
Although a one- dimensional crystal is considered for simplicity, the results in principle apply to a 
two-dimensional lattice if certain conditions are satisfied. The effect of mode-coupling is also briefly 
considered. The relation to previous results on atomic chains, and also to experimental results on 
strongly-coupled dust layers in gas discharge plasmas, is briefly discussed. 
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I. INTRODUCTION 

Dust contaminated plasmas (dusty plasmas, DP) have been attracting significant interest recently. Particularly 
important are dust quasi-lattices, which are typically formed in the sheath region above the negative electrode in 
discharge experiments, horizontally suspended at a levitated equilibrium position at z = zq, where gravity and electric 
(and/or magnetic) forces balance. The linear regime of low-frequency oscillations in DP crystals, in the longitudinal 
(acoustic mode) and transverse (in-plane, shear acoustic mode and vertical, off-plane optical mode) direction(s), is 
now quite well understood. However, the nonlinear behaviour of DP crystals is still mostly unexplored, and has lately 
attracted experimental [1 - 3] and theoretical [1 - 9] interest. 

Recently [5], we considered the coupling between the horizontal (~ x) and vertical (off-plane, ~ z) degrees of 
freedom in a dust mono-layer; a set of nonlinear equations for longitudinal and transverse dust lattice waves (LDLWs, 
TDLWs) was thus rigorously derived [5]. Here, we review the nonlinear dust grain excitations which may occur in a 
DP crystal (here assumed quasi-one-dimensional and infinite, composed from identical grains, of equilibrium charge q 
and mass M, located at x n — nro, n € Af). Ion- wake and ion- neutral interactions (collisions) are omitted, at a first 
step. This study complements recent experimental investigations [1-3] and may hopefully motivate future ones. 

II. TRANSVERSE ENVELOPE STRUCTURES. 

The vertical (off-plane) n— th grain displacement Sz n = z n — zq in a dust crystal obeys the equation |lOllll[ 

+ v + c4,o ( Sz n+ i + Szn-i - 2 5z n ) + uj 2 g Sz n + a {5z n f + [3 (Sz n f = . (1) 

The characteristic frequency 



wt,o = [-qU'(r )/(Mr )] 
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is related to the interaction potential U(r) [e.g. for a Debye-Huckel potential: Uo(r) — (<?/?') e t I Xd , one has 



w o,r> = ^dl exp(-«) (1 + k)/k 3 



(2) 



where u>dl = [q 2 /(-W^c)] is the characteristic dust-lattice frequency scale; Xu is the Debye length; k = tq/Xd is 
the DP lattice parameter]. The gap frequency u g and the nonlinearity coefficients a,j3 are defined via the potential 



pa $(z ) + M 



(3) 



(formally expanded near z , taking into account the electric and/or magnetic field inhomogeneity and charge variations 
|l2|L i.e. leading to an overall vertical force 

F(z) = F el/m (z) - Mg = -d$(z)/dz « -M[to 2 Sz n + a (Sz n ) 2 + (3 (8z n f] + 0[(6z n ) 4 } . 

Recall that F e / m (zo) = Mg. Notice the difference in structure from the usual nonlinear Klein-Gordon equation used 
to describe Id one-dimensional oscillator chains: TDLWs ('phonons') in this chain are stable only in the presence of 
thanks to the field force F e / m (via u> g ). It should be stressed that the validity of this anharmonicity hypothesis is 
indeed suggested real discharge experiments, in particular for low pressure and/or density values, and also confirmed 
by ab initio models 0| (see Fig. [IJ. 

Linear transverse dust-lattice excitations, viz. Sz n ~ cos <fi n (here </>„ = nkr a — tut) obey the optical-like discrete 
dispersion relation [l4| : 



o sin I 



fcro/2) 



2 

: 0J T . 



(4) 



[see that 



The TDLW dispersion curve is depicted in Fig Transverse vibrations propagate as a backward w 
v g ,T — uj' T (k) < 0] - for any form of U(r) - cf. recent experiments [2]. Notice the lower cutoff tOT,min = (u 2 — 4cj^ a) 1 ^ 2 
(at the edge of the Brillouin zone, at k — n/ro), which is absent in the continuum limit, (for k <C r^ 1 ). 
Allowing for a slight departure from the small amplitude (linear) assumption, one obtains: 



Sz n «e(4e*+c.c.) + e 2 



2\A[ 
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Notice the generation of higher phase harmonics due to nonlinearity. The (slowly varying) amplitude = A[e(x 
v g t),e 2 t] obeys a nonlinear Schrddinger equation (NLSE) in the form [7]: 



dA 
dT 



P^ + Q\AfA = 0, 



(6) 



where {X, T} are the slow variables {e(x — v g t), e 2 t}. The dispersion coefficient P is related to the curvature of uj(k) 
as Pt = uj'rf,{k)/2 is negative/positive for low/high values of k. The nonlinearity coefficient 



- 1 f 10a2 
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(7) 



is positive for all known experimental values of the anharmonicity coefficients a, /3 [3]. For long wavelengths [i.e. 
k < k cr , where P(k cr ) = 0], the theory [7] predicts that TDLWs will be modulationally stable, and may propagate in 
the form of dark/grey envelope excitations (hole solitons or voids; see Fig. Qp^b). O n the other hand, for k > k cr , 
modulational instability may lead to the formation of bright (pulse) envelope solitons (see Fig. gp). Analytical 
expressions for these excitations can be found in [7]. 

It may be noted that the modulation of transverse dust grain oscillations clearly appears in numerical simulations 
[H; see e.g Fig. 
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FIG. 1: The (anharmonic) sheath (a) force F(z), and (b) force potential V(z), depicted vs. the vertical distance z from the 
negative electrode, in plasma discarge experiments; figure reprinted from 13]. 

III. LONGITUDINAL ENVELOPE EXCITATIONS. 

The longitudinal dust grain displacements Sx n = x n — nr^ are described by the nonlinear equation of motion Il0| : 

d 2 (5x n ) d(5x n ) 2 

— -^2 I" v — ~T± — = W 0,£ (<^n+l + ox n ^i - 2dx n ) 

-a 20 [(Sx n+ i - Sx n ) 2 - (Sx n - Sx n -i) 2 ] + a 30 [{Sx n+1 - Sx n ) 3 - (Sx n - <5x„_i) 3 ] . (8) 

The resulting linear mode obeys the acoustic dispersion relation: 

uj 2 =4w| )0 sin 2 (fcr /2) = uj\ , (9) 

where lul,o = [U"(ro)/M)] 1 / 2 ; in the Debye case, uj 2 l = 2u> 2 DL exp(-n) (1 + k + k 2 /2)/k 3 . The LDLW dispersion 
curve is depicted in Fig|3J 

The multiple scales (reductive perturbation) technique (cf. above) now yields (~ e) a zero£/i-harmonic mode, 
describing a constant displacement, viz. 

Sx n a e + {u{ 1] e**» + c.c.)] + e 2 {u {2) e 2 ^ + c.c.) + ... . 
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FIG. 2: The TDLW dispersion relation: frequency (square) lj^ vs. wavenumber k. 



The lst-order amplitudes obey the coupled equations [6] : 



1 dT + Fl dX 2 + 1 Ul + 2uo L Ul dX ~ U ' [W) 



d 2 u ( 1] p k 2 d (i) 2 

dX 2 ~ ^-< r^I |Ul 



(11) 



where = w^(fc); {X, T} are slow variables (as above). The description involves the definitions: po = 

—r^U'"(ro)/M = 2020^0 anc l 1o = C^""(' , o)? , o/(2-^) = 3a3o?"Q (both positive quantities of similar order of magni- 
tude for Debye interactions; see in [4, 7]). Eqs. HUH . may be combined into a closed equation, which is identical 
to Eq. ® (for A = u[ 1} , here). Now, here P = P L = Lj'[(k)/2 < 0, while the form of Q > (< 0) prescribes 
stability (instability) at low (high) k. Envelope excitations are now asymmetric, i.e. rarefactive bright or compressive 
dark envelope structures (see Figs.). 

IV. LONGITUDINAL SOLITONS. 

Equation JSJ is identical to the equation of motion in an atomic chain with anharmonic springs, i.e. in the 
celebrated FPU (Fermi- Pasta-Ulam) problem. Inspired by methods of solid state physics, one may opt for a continuum 
description at a first step, viz. 5x n (t) — > u(x,t). This may lead to different nonlinear evolution equations (depending 
on simplifying assumptions), some of which are critically discussed in [9]. What follows is a summary of the lengthy 
analysis carried out therein. 

Keeping lowest order nonlinear and dispersive terms, the continuum variable u obeys |ldj : 

c 2 

u + vix- c L u xx - -^rlu xxxx = -p u x u xx + q iu x ) 2 u xx , (12) 

where (-) x = d(-)/dx; cl = ^L,o r o] Po and 90 were defined above. Assuming near-sonic propagation (i.e. v w cl), 
and defining the relative displacement w = u x , one has 

Wt — aww(; + aw 2 WQ-\- owqqq — (13) 

(for v = 0), where a = po/(2cl) > 0, a — <j<o/(2cl) > 0, and b — C£r 2 /24 > 0. Since the original work of Melands0 
[4], various studies have relied on the Korteweg - deVries (KdV) equation, i.e. Eq. (|13fl for a — 0, in order to gain 
analytical insight in the compressive structures observed in experiments [1]. Indeed, the KdV Eq. possesses negative 
(only, here, since a > 0) supersonic pulse soliton solutions for w, implying a compressive (anti-kink) excitation for 
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FIG. 3: Amplitude modulation of transverse dust lattice oscillations; simulation data provided in the embedded caption; figure 
reprinted from filjjl . 



u; the KdV soliton is thus interpreted as a density variation in the crystal, viz. n(x,t)/no ~ —du/dx = —w. Also, 
the pulse width Lq and height uq satisfy uqLq — est., a feature which is confirmed by experiments [1]. Now, here's a 
crucial point to be made (among others [9]): in a Debye crystal, a w 2a roughly (for k w 1), so the KdV approximation 
(i.e. assuming a ss 0) is not valid. Instead, one may employ the extended KdV Eq. (cKdV) (|13|l . which accounts for 
both compressive and rarefactive lattice excitations (see expressions in [9]; also cf. Fig. 4). 
Alternatively, Eq. I)12[l can be reduced to a Generalized Boussinesq (GBq) Equation 

w ~ vl w xx = h w xxxx + p (w 2 ) xx + q (w 3 ) xx (14) 

(w = u x \ p = —po/2 < 0, q = qo/3 > 0); again, for q ~ go = 0, one recovers a Boussinesq (Bq) equation, e.g. widely 
studied in solid chains. As physically expected, the GBq (Bq) equation yields, like its eKdV (KdV) counterpart, both 
compressive and rarefactive (only compressive) solutions; however, the (supersonic) propagation speed v now does 
not have to be close to cl- A detailed comparative study of (and exact expressions for) all of these soliton excitations 
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FIG. 4: TDL envelope solitons of the (a) dark, (b) grey, and (c) bright type. 
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FIG. 5: The LDLW dispersion relation: frequency ujl vs. wavenumber k (solid curve). We have also depicted: the continuous 
approximation (dashed curve) and the acoustic (tangent) curve at the origin. 



can be found in [9]. 



V. CONCLUSIONS. 



Concluding, we have reviewed recent results on nonlinear excitations (solitary waves) occurring in a (Id) dust 
mono-layer. Modulated envelope TDL and LDL structures occur, due to sheath and coupling nonlinearity. Both 
compressive and rarefactive longitudinal excitations are predicted and may be observed by appropriate experiments. 
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FIG. 6: Bright LDL (asymmetric) envelope solitons: (a) the zeroth (pulse) and first harmonic (kink) amplitudes; (b) the 
resulting asymmetric wavepacket. 




FIG. 7: (a) Grey and (b) dark LDL (asymmetric) modulated wavepackets. 
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relative displacement u x ; (b) grain displacement u. 
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